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9^ Abstract 



We investigate the backflow effect in elementary quantum mechanics - the phenomenon in which 



a- 

qh_ a state consisting entirely of positive momenta may have negative current and the probability flows 



in the opposite direction to the momentum. We discuss various measurement models in which 
backflow may be seen in certain measurable probabilities. We compute the current and flux for 
states consisting of superpositions of gaussian wave packets. These are experimentally realizable 
but the amount of backflow is small. Inspired by the numerical results of Penz et al, we find two 
non-trivial wave functions whose current at any time may be computed analytically and which 
have periods of significant backfiow, in one case with a backwards fiux equal to about 70 percent 
of the maximum possible backfiow, a dimensionless number Cbm ~ 0.04, discovered by Bracken and 



X: 

■ Melloy. This number has the unusual property of being independent of h (and also of all other 

parameters of the model), despite corresponding to an obviously quantum-mechanical effect, and 
we shed some light on this surprising property by considering the classical limit of backflow. 
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I. INTRODUCTION 



A striking but seemingly little-known phenomenon in elementary quantum mechanics is 
the backflow effect. This is the fact that, for a free particle described by a wave function 
centred in x < consisting entirely of positive momenta, the probability of remaining in 
X < may, for certain states, increase with time. That is, the quantum-mechanical current 
at the origin can be negative and the probability flows "backwards" . 

This surprising and clearly non-classical effect, so far unchecked experimentally, was 
flrst noted by Allcock in his seminal work on arrival time in quantum mechanics [l| and 
subsequently highlighted by Bracken and Melloy in 1994, who elucidated some of its features 



In particular, they showed that there is a limit on the total amount of backflow. 
Although backflow means that the probability for remaining in a; < may increase with 
time, the increase can be no greater than an amount Cbm-, a dimensionless number computed 
numerically by Bracken and Melloy to be approximately 0.04. Furthermore, although this 
effect is clearly non-classical, the number Cfym is independent of h (and also of the particle 
mass m and the time duration of backflow). For this reason Bracken and Melloy declared 
Cfem to be a "new quantum number" . 

This remarkable effect has been further investigated by a number of authors. Eveson et al 
signiflcantly reflned the numerical computation of Cfem p]- Similar numerical computations 
were carried out by Penz et al |4] who also gave a rigorous account of the optimization 
problem involved and suggested an analytic form for the state of maximum backflow. Muga 
et al gave some analytic examples of backflow states and explored some aspects relating to 
detection [5]. Berry p\ has also explored various aspects of backflow and related it to the 
phenomonon of superoscillations [3] . Bracken and Melloy have also explored the effect in the 
Dirac equation jsj. Other than these very interesting works, there has been very little work 
on the properties of the backflow phenomenon, although its existence is frequently noted in 
connection with the arrival time problem 

UQ- The purpose of this paper is to explore and 
illustrate various aspects of backflow and in particular to provide concrete analytic examples 
of it. 

We begin in Section 2 with a detailed formulation of the problem. We deflne the current 
and the flux and consider the properties of the spectrum of the flux operator, in terms of 
which the backflow problem is most clearly deflned. In Section 3, we consider some simple 
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measurement models and discuss the ways in which backflow may be seen in the probabihties 
for measurements. In Section 4, we discuss the special properties of the eigenstates of the 
flux operator in measurement models. 

In Section 5 we give some simple examples of states with backflow using superpositions 
of gaussian states. The backflow for such states, is however, rather small. 

In Section 6, we review and repeat the numerical computation of the maximal backflow 
state and eigenvalue. We show that the analytic expression conjectured by Penz et al for 
the maximal backflow state has positive current so cannot be a backflow state. We give two 
non-trivial wave functions with backflow which appear to match very closely the numerical 
solutions for the maximal backflow state by Penz et al {4]. The current at arbitrary times 
of these wave functions is computed analytically and we find that one has a backflow of 
approximately 70 percent of the maximal value. This is a much larger backflow than any 
analytically tractable states previously discovered. 

In Section 7, we consider the naive classical limit — > of backflow, and in particular, we 
address the fact that the bound on backflow Cbm discovered by Bracken and Melloy appears 
to be independent of h. We show that the expected dependency on h reappears in realistic 
measurement models, where measurements are described not by exact projectors but by 
quasiprojectors involving parameters characterizing the imprecision of real measurements. 
Under these conditions the naive classical limit is restored. 

We summarize and conclude in Section 8. 



II. DETAILED FORMULATION OF THE PROBLEM 



A. The Flux 



We consider a free particle with initial wave function iIj{x) centred in x < and consisting 
entirely of positive momenta. Note that it may not of course be perfectly localized in x < 0, 
since this incompatible with positive momenta. We consider the amount of probability flux 
F(ti,t2) crossing the origin during the time interval [^1,^2], defined by 
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where J{t) is the usual quantum-mechanical current at the origin 

,W.«(,.(0.,)f!fi)-fMlM,,o..)) (2.3) 

2m \ ox ox J 

The flux is also easily rewritten in terms of the Wigner function [18] at time t, Wt{p, q), 

F{tu t2)= r dt ! dpdq ^6{q)Wt{p, q) (2.4) 
J m 

It is also useful to write these expressions in an operator form. We introduce a projection 
operator onto the positive x-axis, P = 6{x), and its complement, P = 1 — P = 9{—x). The 
flux may then be written in terms of a flux operator -F(ti, ^2) defined by 

F{UM) = P{h)-P{h) 
dt Pit) 



tl 

t2 



dt ^[H,e{x)] 

I V 



t2 

dt J{t) (2.5) 



tl 



where we have introduced the current operator 



2m 

So Eqs. fl2.2p . (12. 3p may also be written 



J=— {p6{x) + 6{x)p) (2.6) 



F{tut2) = {F{t,,t2)) 

= {Pit,)) - {Pit^)) 
= {Pit^)) - {Pit,)) 



dt {^\Jit)\ij) (2.7) 

tl 

where J(t) = (^|J(t)|^). 

The flux Eq. (12.21) is a difference between two probabilities and is clearly positive when 
those probabilities behave according to classical intuition, i.e., when the probability of re- 
maining in a; < decreases monotonically. For this reason, the flux is often proposed as 
the provisional semiclassical answer to the arrival time problem: what is the probability for 
crossing the origin during the time interval [ti,t2]? This is discussed at length elsewhere 



4l6l | and, although the arrival time problem forms the backdrop to the current work. 



the properties of the flux pose an interesting set of problems in themselves and it is these 
problems we focus on. 

As indicated already, in the full quantum-mechanical case, the flux can be negative. The 
above formulae give some clues as to why this is the case. First of all, since the Wigner 



function can be negative [18|], Eq. (l2.4p suggests that the flux can be negative for certain 
states. More precisely, negative Wigner function is a necessary condition for negative flux, 
which indicates that it relates to states with interferences in position or momentum. Negative 
Wigner function is not a sufficient condition since the integral in Eq. ( 12.41) may yield a positive 
expression, even for negative W. 

The second clue to the possible negativity comes from the current operator Eq. fl2.6p : the 
two operators p and 6{x) are non- negative (on states with positive momentum), but since 
they do not commute, the current operator J is not a positive operator. 

B. Most Negative Flux as an Eigenvalue Problem 



Following Bracken and Melloy [2^, a useful way to find the states with negative fiux is 
to look at the spectrum of the fiux operator Eq. (12.51) . We thus look for the solution to the 
eigenvalue problem 

F(ti,t2)|<f) = A|$) (2.8) 

where the states |$) consist only of positive momenta. (We choose an opposite sign conven- 
tion to Bracken and Melloy which means that the backfiow states have A < 0). The most 
negative value of the fiux F{ti, ^2) is then the most negative eigenvalue of the fiux operator. 

It is convenient to choose the time interval [ti, ^2] to be [— T/2, T/2], as is easily achieved 
by time evolving the state, and the eigenvalue equation in momentum space then reads 

1 f°° „ sin[(p2 - fc2)T/4m^] , 

- dk — f-^ ^ ^k) = X^{p) (2.9) 



We then define rescaled variables u and v hy p = 2^Jmh/Tu and k = l^fmKjTv and the 
eigenvalue equation is then 

1 ('^ sin(w^- t;^) ^^^^ ^ ^^^^^ ^2.10) 

TT Jo (« - 

where </)(«) = (m^/4T)^/^$(p) and is dimensionless. Note that all physical constants /i, m, T 
have dropped out of this equation so that the eigenvalues A are dimensionless and indepen- 
dent of m and T. Note also that the equation is real so one may take the eigenstates 



(f){u) to be real. It is also useful to record the result that the flux for any state (f){u) in these 
variables is given by 

F(-T/2,r/2) = - / du dv(t)*{u) -^ r^(t>{v) (2.11) 

Jo Jo [u - v) 

The eigenvalue equation Eq. fl2.10p clearly has approximate solutions with eigenvalues 

close to 1 or consisting of wave packets which cross the origin either well inside or well 

outside the interval [— T/2,T/2]. Further study of this eigenvalue equation has been carried 

out by a number of authors, both numerically and analytically 2|-|^. The eigenvalues lie in 

the range 

-Cfe^<A<l (2.12) 
where ctm was computed numerically and found to be 

Cft^^ 0.038452 (2.13) 

It was conjectured in Ref. 21] that the spectrum is discrete in the interval [— c;,m,0] but 
continuous in the interval [0, 1]. Penz et al 4|] asserted, from their numerical studies, that 
the extremizing state has the form sinw^/w. However, as we shall see, this can be true at 
best asymptotically, since the current for this state can in fact be computed analytically and 
is positive. We will, however, flnd good analytic expressions approximating the numerical 
results in what follows. 

At present there is no analytic account of the properties of the results Eqs. fl2.12p . (12.131) . 
Physically, backflow is related to the fact that specifying both positive momenta and asking 
for the probability of remaining in x > correspond to incompatible measurements in 
quantum mechanics. That is, it is related to the fact that the operators 6{p) and 6{—x) 
do not commute. This leads to the question, is there an analytic calculation of the most 
negative eigenvalue —Cbm, perhaps involving the non-commutativity of 6{p) and 9{—x)l We 
do not have an answer to this question but it remains an important question for future 
study. 

The fact that the eigenvalues are independent of h means there are some potential prob- 
lems with the naive classical limit and we address this below. The eigenvalues are also 
independent of T. This simply means that the duration of a period of negative current can 
be arbitrarily long, as long as the total flux over that time period is bounded from below by 



III. BACKFLOW AND MEASUREMENT MODELS 



In this section we relate the above results on backflow to measurements. This is partly 
to begin to address the practical question of how backflow is measured, but also to get some 
insight into the negativity of the flux. 



A. Explicit Measurement of Backflow 

The flux Eq. fi2.2p is defined as the difference between two probabilities, therefore the flux 
can be measured by measuring these two probabilities. This requires two ensembles, each 
prepared in the initial state On one ensemble, measurements are made to determine if 
the particle is in a; > at time ti, hence determine (P(ti)). On the second ensemble, the 
same measurements are performed at time t2, which thus determines {P{t2))- From the two 
results the flux can determined. (Further discussion of the possible measurement of backflow 
is given by Bracken and Melloy [2|). 

Note that the above is not the same as performing sequential measurements of position 
on the same ensemble. We will come to these sorts of measurements below. 



B. A Simple Measurement Model for Arrival Time 

Some mathematical insight into the negativity of the flux may be obtained by looking 
at simple models for measuring the arrival time. This is because such models, if properly 
constructed, yield a positive probability which will, however, be approximately the same as 
the flux in some limit. Hence by comparing such models with the flux we may be able to 
see the origin of the negativity. 

The simplest model for measurement of the arrival time involves simply measuring to see 
if the particle is in a; < at time ti and then in a; > at time ^2- This probability is given 
by 

p{h,t2) = {^P\Pih)Pit2)PihM (3.1) 

which is clearly positive. It gives a simple notion of arriving at the origin during the time in- 
terval [^1,^2] (but ignoring issues about multiple crossings). Using the flux operator Eq. fl2.5p 
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this may be rewritten 

= dt {ilj\P{h)J{t)P{ti)\ij) (3.2) 

This coincides with the flux Eq.f l2.7p except for the projection operators onto x < at ti. 
Since Eq. (13.21) is positive, this means that the negativity of the flux comes entirely from the 
part of the state which is aheady in x > at the initial time ti. 



C. A Complex Potential Model for Arrival Time 



involving; model detectors, for 
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19l-|22l|. We again consider 



More elaborate approaches to the arrival time problem 
example) naturally lead to the use of a complex potential 
an initial wave packet starting in x < with positive momentum and seek the arrival 
time probability distribution n(r)(ir for crossing the origin between r and r + dr. We 
consider a complex absorbing potential of step function form in x > so the Hamiltonian 
is i^o ~ '^Vo^(^'), where Hq is the free Hamiltonian. We deflne the survival probability N{t) 
to be the norm of the state at time r after evolution with this complex Hamiltonian. The 
arrival time distribution is then given by 

dN 



H r 



dr 

2Vo(V^|e(*^°-^°^("»"^(x)e(-*^°-^°^(^))"|V^) (3.3) 



We seek a simple form for this expression which exposes its dependence on the current 
operator and thus gives some idea as to how it will be affected when backflow is present. 
Differentiating with respect to r, we get 

^ = -2K)H + 2\/o(V'|e(*^o-^°^(^))Ve(-*^o-^«^(^))"|^/^) (3.4) 
dr 

where J is the current operator Eq. (12.61) . Eq. (l3.4p is a differential equation for H(r) which 
is easily solved to yield 

H(r) = 2Vo r dt e-2^»(^-*) {^l;\e('^o-Voeimj^{-iHo-Voemt\^) (3.5) 
J —00 

where we have assumed that H(r) — )■ as r — t- —00. Eq. (l3.5p is the exact expression for H(r) 
and displays the dependence on the current operator J. It is positive by construction, even 

8 



though j is not a positive operator. The probabihty for crossing during the time interval 
[ti, ^2] then is 



t2 



dt n(t) 



(3.6) 



This is the analogue of Eq. fl3.2p . 

It is not easy to see how the presence of backfiow states in Eq. fl3.5p may register in 
the probability distribution n(r). The expression is, however, simpler in the usual weak 
measurement approximation (small Vq), which involves neglecting the complex potential 
terms in the bracket expression, yielding 



(3.7) 



where \ipt) = e ''^"^lip). This is the expected semiclassical result |23|, |2| (although 

note that the new derivation given here is considerably shorter than those given elsewhere!). 

The quantity n(r) corresponds to the arrival time distribution measured by a realistic 
measurement so we in principle be determined ex per imentally. The current can in principle 
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2^. 



be extracted from Eq. fl3.7p by deconvolution 

Note that Eq. fl3.7p is not necessarily positive, due to the negativity of the current in 
certain states. The positivity has been lost in going from Eq. fl3.5p to Eq. fl3.7p because we 
took the limit Vq — ?■ in the bracket expression but not in the rest of the expression. 



IV. MEASUREMENT PROBABILITIES FOR FLUX OPERATOR EIGEN- 
STATES 

Simplifications to simple probabilities of the form Eq.( ]3.ip occur in the particular case 
in which the initial state is an eigenstate of the flux operator. Since such states include the 
state with maximal backfiow, these states are particularly interesting to study. 

We write the flux operator eigenvalue equation Eq. fl2.8p as 

(P2-Pi)|$a) = A|<I>a) (4.1) 

where we have introduced a simplifying notation P(t2) = P2, P(ti) = Pi. The eigenstates 
|$a) consist of entirely positive momenta. They are localized in a; < in the distant past 
and in X > in the distant future j^. In the following, we are most interested in the 
eigenstates with negative eigenvalues, for which, as mentioned, the spectrum is conjectured 
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to be discrete 2| and therefore the norm ($a|$a) exists. (Indeed, we will see later that the 
conjectured asymptotic form of the state with most negative eigenvalue is normalizable). 
For eigenstates belonging to the continuous part of of the spectrum (which are probably 
the states with positive eigenvalues, which are of less interest) it is necessary to construct 
normalizable states of the form 

m = j dX c(A) 1$^) (4.2) 

for some smearing function c(A) peaked around a particular A. We will assume that this has 
been done where necessary in what follows. It will not effect our essential conclusions which 
in any event mainly concern the allegedly discrete negative eigenvalues. 
Averaging Eq.f l4.ip yields 

(P2) - {Pi) = A (4.3) 
Multiplying by P2 and averaging, and by Pi and averaging yields 

(P1P2) = (1-A)(P2) (4.4) 
(P2P1) = (l + A)(Pi) (4.5) 

Since the right hand side of both of these equations is real we have that (P1P2) = (P2-P1) 
and combining with Eq. fl4.3p we find 

(^1) = (A) = ^(l-A) (4.6) 
{P2) = (A) = ^(l + A) (4.7) 

Proceeding in this way we can easily find the averages of any string of projectors in a 
flux operator eigenstate. In particular, we find the probabilities 

(P1P2P1) = i(l + A)(l - A)^ = 1(1 - A)(l - A^) (4.8) 

(AP2P1) = 1(1 - A)(l + A)^ = ^(1 + A)(l - A^) (4.9) 

(P1P2P1) = ^(1-A)A2 (4.10) 

(P1P2P1) = ^(1 + A)A2 (4.11) 

These four probabilities represent the two ways of crossing the origin and the two ways of not 
crossing the origin during the time interval [ti, t2], as characterized by position measurements 
onto X > and x < at the two times. 
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It is also useful to compute an expression which represents the interference between 
crossing and not crossing, 

{P,P,P,) = \\{l-\') (4.12) 

There are three more similar terms but they have similar characteristics. This interference 
ierm is best understood in terms of the decoherent histories approach to quantum theory 
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25| but it is not necessary to go into that here. Basically, the smallness of this term is 



a measure of the degree to which classical intuition can be relied on. This interference term 
is zero when A = or A = 1, so only in these cases can we expect classical intuition to be 
reliable. Hence backflow states with A < are ruled out, as expected, but note that positive 
values of A between and 1 may also be non-classical, by this criterion. 

Consider now the probabilities Eqs.f l4.8p -f l4.lip . The two crossing probabilities Eqs. (14.101) . 
(14. lip have the property that 

(P1P2P1) 1-A , , 

- = r 4.13 

which is greater than 1 when A < 0. That is, for backflow states, the probability of crossing 
from right to left, is greater than the probability for crossing from left to right. Recalling that 
the overall motion of the wave packet is from left to right, the backflow states do indeed have 
the non-classical property of flowing backwards, here deflned in terms of sequential position 
measurements. 

Another simple observation from Eq. fl4.8p is that 

i^^ = (l-A)^ (4.14) 

which is greater than 1 when A < 0. This means that the probability of flnding the particle 
in X > at time ^2 is increased if an extra measurement is inserted at time ti, checking to see 
if the particle is in x > 0. This is surprising since one would expect the extra measurement 
at ti to restrict the flux of wave function entering x > at t2, yet it increases it. This 
suggests that there is some sort of diffraction effect involved - the extra projector at ti 
diverts extra flux into a; > at ^2- 
Finally, note also that 

= {l + Xf (4.15) 



which is greater than 1 for < A < 1. (We exclude A = 1 since both numerator and 
denominator vanish.) This is surprising for the same reason the previous case is and also 
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suggests some sort of difFraction effect. Although note this happens for positive A, confirming 
the earher statement that values of A in this range may also give non-classical states. 



V. BACKFLOW FOR SUPERPOSITIONS OF GAUSSIANS 

n 

Bracken and Melloy gave two explicit examples of states displaying backfiow [2] . Although 
these examples served to demonstrate the existence of the effect, the particular states chosen 
were rather unphysical. In this section we show that the backfiow effect also arises in the 
more familiar, and also potentially experimentally realisable, setting of a superposition of 
two gaussian wavepackets. A single gaussian has positive Wigner function must have positive 
current. Eg. (12.41) . The Wigner function of a superposition of Gaussians may, however, be 
negative in some regions, so the current can be negative. 

Of course the problem with using gaussian wavepackets is that they have support on 
both positive and negative momentum, and so as well as demonstrating the appearance of 
negative current, we will also have to show that this is not the result of any initially negative 
momentum. We will see that whilst superpositions of gausian states do indeed give rise to 
backfiow, the size of the effect is considerably smaller than the theoretical maximum. 

We begin with the simple case of a superposition of two plane waves, as considered in 



.0. 



Ref.l2l This can be turned into a more physical state by replacing the plane waves with 
gaussians tightly peaked in momentum, without affecting the basic conclusion that the state 
displays backfiow for well chosen values of the various parameters. In this Section we work 
in units in which h = 1 and we set the particle mass m = 1. 
We start with the unnormalised state. 



i)'{x,t) = ^ Akexp[ipk{x - pkt)] (5.1) 

fc=l,2 



where the Ak are taken to be real. (One could of course add to each component an arbitrary 
phase, but this is an unneccesssary complication.) The current at the origin for this state is 
given by, 

J{t) = Ajpi + AIp2 + AiA2ipi + P2) cos[(Ei - E2)t] (5.2) 

This oscillates between a maximum value of {Aipi + A2P2){Ai + A2) and a minimum value of 
{Aipi — A2P2){Ai — A2). Thus, for instance, if Ai > A2 and Aipi < A2P2 this state displays 
backfiow. 



12 



Consider now the normalised state 

= E A,^=i===exp (^p,ix-p,t) - ^{^^) . (5.3) 

This is a sum of two initial gaussian wavepackets with equal spatial width a, evolved for a 
time t. If we let a — ?■ oo we essentially recover Eq.f l5.ip . The idea is that if we take a to 
be large enough, the current at the origin is the product of Eq. (15.21) and a slowly varying 
function, so that the conclusions about backflow still hold. The analytic expression for the 
current is somewhat long and complex, and we will not give it here. Instead we show below 
two sets of plots of the current at the origin and the probability of remaining in x < 
as functions of time for the state in Eq. fl5.3p for the following two sets of parameters. In 
Figs. ([ID and ([2]) we plot the parameters 

p, = 0.5, P2 = 2, a = 10, Ai = 1.7, A2 = 1. (5.4) 

and in Figs. ([3]) and (jlj) we plot the parameters 

pi = 0.3, P2 = 1.4, a = 10, A, = 1.8, A^ = 1. (5.5) 

We clearly see from these plots that there are several intervals during which the current is 
negative. These examples show that the backflow can occur in several disjoint time intervals. 
A magnification of one of these backflow region is shown in Fig.(l5l). The effect is robust with 
respect to small changes of the parameters, Eq. (l5.4p or Eq. (l5.5p . which were chosen because 
they give reasonable amounts of backflow. 

The set of parameters in Eq. (15.50 give rise to the greatest amount of backflow we have 
been able to find, although we have not performed anything like a comprehensive search of 
the parameter space. The value of the flux during the largest period of backflow is 

F = [ dtJ{t), (5.6) 
Jti 

where J(t) is the current, Eq. (l2.3p . and the interval [^1,^2] is chosen such that the current is 
negative for the whole of this time. Computing the flux during the largest backflow interval 
gives, 

F ^ -0.0061, (5.7) 
or about 16 percent of the theoretical maximum. 
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FIG. 1: Plot of the current for a wave- 
function consisting of a superposition of 
two gaussians, with the parameters given in 
Eq.dSai). 
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FIG. 3: Plot of the current for a wave- 
function consisting of a superposition of 
two gaussians, with the parameters given in 
Eg. ([53]) . 



FIG. 2: Plot of the probability for remaining 
in X < for a wavefunction consisting of 
a superposition of two gaussians, with the 
parameters given in Eq. (j5. 




FIG. 4: Plot of the probability for remaining 
in X < for a wavefunction consisting of 
a superposition of two gaussians, with the 
parameters given in Eq. ()5.5p . 



It is important to check that this probability backflow cannot be explained by the tiny 
probability of having negative momentum which comes from this gaussian state. An order 
of magnitude estimate will suffice here. We have two gaussian wavepackets centered about 
different momenta. Consider the wavepacket centered around p = 0.3. The probability 
that a measurement of the momentum of this state would yield a negative answer is given 
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FIG. 5: Close up of Fig HI The probability 



is clearly seen to increase between t ~ 2 and 



t 



4. 



approximately by, 




~ 10 



10 



(5.8) 



J —CO 



The negative flow of probability is therefore entirely due to the backflow effect. 

VI. AN APPROXIMATION TO THE BACKFLOW MAXIMISING STATE 

Backflow states may be found by solving the eigenvalue equation Eq. fl2.10p . The numer- 

I 

ical work of Penz et al ^] yielded a plot of the approximate eigenstate satisfying Eq. fl2.1Up 
giving the most negative eigenvalue —Cbm, i-e., the largest amount of backflow. They also 
conjectured the approximate analytic form of this state. However, we have found that the 
current of this state can be computed analytically and is positive so cannot be the optimizing 
state, except asymptotically. It is therefore still a problem of great interest to find analytic 
expressions for wave functions with backflow which match the numerical results as closely 
as possible and for which the current at arbitrary time may be computed analytically. This 
is what we do in this section. We will not find analytic solutions to the eigenvalue problem, 
Eq.f l2.10|) . Rather, inspired by the numerical solutions to the eigenvalue problem, we will 
exhibit analytic expressions for wave functions, compute their current at time t analytically, 
and show that they have significant negative flux, calculated using Eq. fl2.1ip . 
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A. Numerical results 



We first review the numerical results of the computation of the backflow state. We have 
repeated the numerical analysis of Penz et al. {4] of the optimizing state and its current. This 
is purely for comparison with our analytic results and we do not claim any improvements in 
accuracy over Penz et al. 

Our numerically computed maximum backflow state has the asymptotic form, 



4>as{u) 



sm -u 



cos M 



(6.1) 



u u 

In Fig. ([6]) we plot the numerically computed maximum backflow state together with (pas for 
a = 0, b = —0.1. Comparing by eye it seems like these parameters produce the best fit, but 
we will consider states with the more general form Eq.( ]6.ip . We plot in Fig.Q the current 




FIG. 6: Plot of the backflow maximising state (solid line) together with (pas for a = 0, 6 = — 0.1 
(dashed line). 



computed from the numerically obtained backfiow maximizing state. These two plots are in 
agreement, in general shape, with the numerical results of Penz et al J] and we will compare 
our analytic results with these plots in what follows. 

Penz et al ^\ suggested for the approximate form of the backfiow state the function 
0(m) ~ sin(-u^)/M. This cannot be correct, since it results in a positive current at t = 0, 
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FIG. 7: The current, J{t), as computed from the numerically obtained backflow maximizing state. 

contrary to the numerical result shown in Fig. ([7]). We see this as follows. In terms of the 
dimensionless quantities introduced in Section 2, the current at t = is 

1 f°° 

J(0) = — / dudv{u + v)(j){u)(j){v) 
2vr Jo 

2\ /"CXD 







U 



sin[u J I T . , 2\ 
au / dvsin[v ) 



1 

-a/- > 0. 



(6.2) 



However this does not rule out the possibility that the backflow maximising state has the 
form smu^/u for large u (even though Fig. ([6]) would seem to suggest otherwise). 

B. Guess 1 



Our first guess is the momentum space wavefunction (piiu) given by, 
0i(n) =iV 



{\-C{u)) + a{\-S{u)) 



a e 



where N is some normalisation factor. Here 



C{u) = FresnelC | y —u 



dx cos (a; 



(6.3) 



(6.4) 
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and similarly for S{u). This state has the asymptotic form 
Note that, 



sinfw^) cos(u^) 
— ^ — - + a- 



u 



u 



dzu {cos^z^v^) + asin(z^u^)) 



which is a form we shall use below. Note that the norm is given by, 

= -1=(1 + + 2a(l + a){V2- 1)) 



(6.5) 



(6.6) 



(6.7) 



We plot the wavefunction in Fig.([H]), and we see that it shows reasonable agreement with 
the numerical result. 




FIG. 8: (i){u), for a=0.4 (solid line), with the exact numerical result for comparison (dashed line). 



We wish to evaluate the flux, Eq.( l2.11| ) which may be written 



dt J{t) 



where we have introduced the current 



J{t) = — / dudv{u + v) eyiY){it{u —v ))(f){u)(f){v) 
2vr Jo 



(6.8) 



(6.9) 



written in terms of the dimensionless variables introduced in Section 2 (and here t is a 
dimensionless time parameter). 
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Because we are working only with approximate eigenstate of the flux, the interval during 
which the J{t) is negative may not coincide exactly with [—1,1], so for that reason, we take 
the flux instead to be 

F= dt J{t) (6.10) 

and adjust the interval [^1,^2] to match the region of negative flux. Eg. (16. 9 p can also be 
written as, 



TT 



J{t) = Re I — / duexp{itu'^)(l){u) / dv v exp{~itv'^)(f){v) 



TT 



-Re{U{t)V{t)) 



(6.11) 



so that we can compute the current by first computing each of the integrals in Eg. (16. lip 
seperately. 

We begin by computing the U integral in Eg. (l6.1ip . We use Eg. (16.60 to write 



/2 pcfo poo 
- du dz- ((1 

(6.12) 

We would like to change the order of integration at this point, but we cannot, since the 
u integral only converges conditionally. To remedy this we introduce a convergence factor 
exp(— e-u^) where e > 0. We can then write. 



00 POO 



X 



(1 — ia) exp{i{t + z )u — eu ) + {1 + ia) exp{i{t — z )u — eu ) 



N 



2\/27r Ji 

iN 



dz 



1 — ia 



1 + ia 



2V27r 



[1 — ia) 



e-i{t + z'^) e-Ot-z"^ 
ArcTan(;7^) ArcTan( — 



N 



2J2Ti{t + ie) 



Tx{l + i){l + a) 



- (1 + ia) 



[i + a)ArcTan( 



a/— t — ie 
1 



z=l 



y/t + ie' 



— {i — a)ArcTanh(- 



'^Jt + if:' 

where e{t) is the signum function, and we have used the standard integrals, 

ArcTan(^) 



/du 
a + 



du- 



ArcTanh(^) 



a — u'^ 



a 



(6.13) 

(6.14) 
(6.15) 
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We now turn to the V integral in Eq.f l6.11|) . Again, we use Eg. (16. 6 p to write, 



oo ^2 



X 



dv / dz — 
Ji 2 

((1 - ia) exp{-i{t - z^)v'^) + (1 + ia) exp{-i{t + z'^)v'^) 



(6.16) 



As before, in order to change the order of integration we insert a convergence factor. 



N 



(1 — ia) exp{—i(t — z^)t>^ — et>^) + (1 + ia) exp{—i(t + ^ — ef ^) 
1 — ia 1 -\- ia 



4^2 ii 



4^2 

AT 



(e + Z(t - 22))3/2 (e + ,(t + ^2))3/2^ 

z^i\/ + t — ie 



(1 — 'ia)^- . , , „ : — (1 + za 



(t-ze)(22 



(t-ze)(22 + t_ie) 



J 1 



''1 — ia) 



t — ie 



1 - 



1 



'1 + m) 



t — it 



1 



(6.17) 



where we have used the standard integrals, 

1 



du 



u 



1 u 



(6.18) 
(6.19) 



(a -1x2)3/2 

Given U and we may now plot the current, Eq. (l6.1ip . The plot is shown in Fig.(|9]). 
It has a substantial period of backflow and is in broad agreement with the numerical result 
Fig. ([7]), although differs significantly in the behavior near t = ±1. We have chosen the 
value a = 0.4 which approximately maximises the backflow for this wavefunction. The flux, 
Eq. (l6.10p . for this choice may be calculated by numerically integrating the current and is 
approximately 

F = -0.02095 (6.20) 

The amount of negative flux obtained is of the order of 55 percent of Cbm, which is a much 
greater fraction than we were able to achieve in Section 5 using superpositions of gaussians. 
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FIG. 9: Current, Ji(t), for a = 0.4 and e = 10"^. 



C. Guess 2 



Our second guess is the momentum space wavefunction 



ae-"" + (- - C{u)) 



a,b & 



where is a normahsation factor. This has the asymptotic form, 

,sin(M^) 



u) ~ N- 



u 



(6.21) 



(6.22) 



We plot 02 in Fig.ffTOl) for the values of a and h which maximize backflow. We see good 
agreement with the numerical result. 
We first calculate the norm N , 



N- 



du ( a^e-^^" + (- - C{u)f + 2ae-''"(- - C{u)) 



(6.23) 



The first two terms can be computed easily enough. The last term is more challanging. 
Using the standard integral. 



COS 



2 Vvr 



sm 



(6.24) 
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FIG. 10: for a=0.6 and b=2.8 (solid line), with the numerical result for comparison (dashed 

line). 



we find, 



= — + 



a 

+ T 



2b 40F b 



2a 
T 



TX 



COS I — I - 



sm,- 



(6.25) 



We now turn to computing the U integral in Eq.( l6.11|) . We wish to compute 



N 
aN 



due^ 



ae 



-bu 



+ - C{u] 



due-^'-''^''\^-C{u)) 



(6.26) 



2 

where we have added a convergence factor e~^^ . The second integral we have already seen, 
it is just Ui{t) for out first guess of 4>i{u), with the coefficient of the term involving S{u) 
taken to be zero. The first integral can also be done easily, using the standard integral 



1 /TT /S"^ 

dx exp{—ax'^ — (3x) — 2y ~c4^Erfc 



/3 



2v^ 



for Re(a) > 0. 



(6.27) 



Combining these gives. 



N 



2 V ^ + ^e 



2^2Ti{t^ie) 



7r(l + i) 



iArcTan( 



zArcTanh( 



(6.28) 
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Next we compute the V integral in Eq.f l6.11|) . 

POO / 1 

V2{t) = N dvve-^'+''^''' iae-'"' + -- C{v) 



oo 



= aN I dwe-'^-^'^''^"' + [ dvve-^'+''>\l-Civ)) (6.29) 
Jo Jo 2 

As with U2{t), the second integral can be simply written down by comparing to Vi{t) for our 

first guess of (pi{u). The first integral is also easily performed, using the standard integral, 

„ <^-e-'-* = ^(l-fy^exp(-gErfc(^)). to. Re(„) > 0. (6.30) 

We thus find, 



V2it) = ( 1 - ^ J^exp ( -— -- ) Erfc 



2{e + it)\ 2\l e + it \ 4:{e + it) J \2^7Tit 



N 

4 



4\/2 



i ( ^ 1 \ Vi 1 



(6.31) 



t — ie \ Vl — t + / t — ie\ v^l + t — ie. 

Now we have U and V we can compute the current Eq. 06.111) and, by numerical integra- 
tion, the fiux Eq. fl6.10l) . The maximum amount of negative fiux we can generate is 

F = -0.02757, (6.32) 

which occurs for the parameters a = 0.6, b = 2.8. This corresponds to about 70 percent of 
the maximum Cf,m- We plot the current J{t) for these parameters in Fig. (ITT]) . The current is 
closer to the numerical result. Fig. ([7]) than our first guess. Fig. ([2]) but still lacks the correct 
behavior as t — ?■ ±1. 

Finally, to give a clear illustration of the backfiow phenomenon, in Fig. (fT2l) we plot the 
probability that the state will be found in x < as a function of time, for the state Eq. fl6.2ip . 
The probability decreases over the whole interval plotted, but has a very noticeable period 
of increase between t = ±1. 

VII. THE CLASSICAL LIMIT OF BACKFLOW 

Some insights into the properties of backfiow may be found by looking at its classical 
limit. In the usual account of emergent classicality for the free particle considered here, 
one considers a larger system in which the particle is coupled to a wider environment and 
one considers the evolution of the reduced density matrix of the particle only Q . It is 
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FIG. 11: Current, J2{t) for a = 0.6, b = 2.8, and e = 10"'^. 



well-known that, in a wide variety of such open system models, the Wigner function of the 
particle will become positive after a short period of time 28|. From Eq. fl2.4p . it is then easily 
seen that the flux will then be positive, so the backfiow clearly goes away in the standard 
approach to the classical limit. A detailed discussion of the arrival time problem in the 
presence of an environment was given in Ref. 29| and in this work the resulting positivity of 
the current, after finite time, is clearly seen. 

However, there is a more interesting and subtle question, noted by Bracken and Melloy 
{2!, which is that the eigenvalues of the flux operator are independent of h, as we saw in 
Eq. (12.1 op . despite the fact that the existence of negative eigenvalues (negative flux) is clearly 
a quantum phenomenon. This means that in the naive classical limit, 0, the backfiow 
does not go away, as one might expect. Of course, this "limit" is an oversimplification of 
what the classical limit means, but despite this, it is still disconcerting that this obviously 
quantum phenomenon is apparently independent of h. 

This situation is reminiscent of another situation without naive classical limit, namely, 
scattering off a step potential, where it is known that the quantum-mechanical reflection 
coefficient is independent of h, so does not go to zero as /i — ?■ 0, contrary to classical 
expectations. The origin of the problem is the use of an exact step potential. If instead a 
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FIG. 12: Probability P{t) that the state will be found in x < at time t for the wavefunction (p2- 
The probability is clearly seen to increase during the interval [—1, 1]. 



smoothed off step is used, with lengthscale a describing the size of the smoothirig region, 
then reflection does indeed go away if /i — > with a held constant and non-zero 30|]. The 
point here is that the exact step potential is an idealization that fails to capture all physical 
properties. Replacement with a more realistic potential restores the naive classical limit. 

In the backflow situation, we may therefore also expect to get a reasonable naive classi- 
cal limit by small modification of the situation. In particular, instead of defining the flux 
operator in terms of exact projection operators P = 6{x), we define it in terms of a quasipro- 
jector Q. This seems reasonable since, as discussed earlier, backflow can be measured by 
measuring whether the particle is in x > at two different times and, due to the inevitable 
imprecision of real measurements, such measurements are best modeled by quasiprojectors. 
A convenient choice of quasiprojector is 



Q 



dy 6a{x - y) 



(7.1) 



where 5a{x — y) is a smoothed out 5-function, 

1 



5^{x - y) 



(27ra2)i/2 



exp 



- yf 

2a2 



(7.2) 
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This goes to the usual 5-function as a — )■ and then Q — )■ 9{x). If we replace P with the 
quasiprojector Q in the expression for the flux, we get Eg. (12.71) but with the current operator 
replaced by 

J = ^{p5„{x) + 5^{x)p) (7.3) 

The resulting flux will, loosely speaking, by less negative, since the commutator between p 
and 5fj{x) becomes smaller as a becomes larger. 

With the quasiprojector, the flux written in the form Eg. (12. 11]) acquires an exponential 
factor 

F(-T/2,T/2) = - / du dv (tf{u) ""^ ) ^~aHu-vf ^^^^ ^^ ^^ 

where the dimensionless number a is given by = 2ma'^/hT and the eigenvalue equation 
Eq. (l2.10 l ) will acquire the same exponential factor, 

1 f°° , sin(M^ — f^) „2/,, ,,-,2 ,/ X , ^ 

dv — ^ e"" (f){v) = A0(m) (7.5) 

[u — v) 



TT 







This means that the eigenvalues A will now depend on a, so we write A = A (a). Through a 
they will therefore depend on h and the "limit" ^ — ?■ now clearly means the regime a ^ 1, 
that is, h ^ Ima"^ jT . Hence, in a more realistic measurement situation, the bound on the 
total backflow - the most negative eigenvalue of Eq. (l7.5p - will depend on U and the limit 
— )■ may now be more meaningful. (Bracken and Melloy noticed a similar phenomenon 
in the context of the Dirac equation, where the presence of the speed of light as another 
physical parameter permits the construction of a dimensionless parameter analogous to a 
above [s'].) 

A reasonable conjecture is that the negative eigenvalues will increase with a and also that 

A(a) > -Cfom (7.6) 

for all a, so that the Bracken- Melloy bound — emerges as a lower bound on the eigen- 
values, achievable only in the limit a — > 0. It seems unlikely, however, that all the negative 
eigenvalues will all become positive or zero, except perhaps in the limit a — )■ oo. This 
behaviour is best explored numerically, which we now consider. 

We consider the behaviour of the most negative eigenvalue A (a) of Eq. (l7.5p . We have 
not been able to solve this equation analytically, so instead we have obtained numerical 
estimates for A(a) for various values of a, and we plot the result in Fig.( IT3|) . The value of 
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A (a) does indeed increase with a, tending to zero asymptotically. In fact, numerical solutions 
are consistent with the asjTiiptotic form, 



A(a) 



(7.7) 



for large a. We have not managed to prove this analytically however. 




FIG. 13: Plot of the size of the most negative eigenvalue of Eq. ()7.5p . A(a) as a function of a. 



To confirm that this result about the negative eigenvalues is significant, we need also 
to compare it with the behaviour of the positive eigenvalues for large a. (It could be, for 
example, that they all go to zero, rendering the above result spurious). As noted in Section 
2, wave packets which clearly cross or do not cross the origin during the interval [— T/2, T/2] 
are approximate eigenstates of Eg. (12.101) with eigenvalues 1 or 0. It is reasonably clear that 
they will also be approximate eigenstates of the modified eigenvalue equation Eq. (l7.5p as 
long as the wave packets cross the origin sufficiently far from the end-points of the interval 
[—T/2, T/2], since under these conditions the incoming wave packet does not notice the 
smearing of the projector into a quasiprojector. This is also backed up by numerical work. 
We have computed a number of eigenvalues for different values of a. There appears to be 
a reasonably even distribution of positive eigenvalues in the interval [0, 1] for a range of 
values of a, so, unlike the negative eigenvalues, this part of the spectrum is not significantly 
changed as a becomes large. 
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In summary, for smeared projectors, Eq. fl7.ip . which are perhaps better models for real 
measurements than exact projectors, the eigenvalues and in particular the lowest eigenvalue 
representing most negative flux are dependent on h (and m and T). Numerical work indicates 
that the most negative eigenvalue increases with a. It is still negative for finite a, indicating 
that backflow will still be present for more realistic measurements. The lowest eigenvalue 
appears to go to zero for a — )■ oo. This indicates that all the negative eigenvalues go to 
zero (or become positive) in the naive classical limit /I — )■ 0. (This is in contrast to the 
more realistic classical limit of a particle coupled to an environment, where it can be shown 



that the backflow effect does indeed vanish exactly after a finite time 
positive eigenvalues are not significantly affected. 



29|). By contrast the 



VIII. SUMMARY AND CONCLUSIONS 



The purpose of this paper was to explore and illustrate the backflow effect in a variety 
of different ways. After setting up the problem in Section 2, in Section 3 we discussed 
measurement models that exposed certain aspects of backflow. Eq.(3.2) establishes that 
backflow arises from the part of the state which is already in x > 0. The complex potential 
model of Section 3(C) corresponds to a number of reasonable realistic measurement models. 
The current can be obtained from the measured probability by deconvolution, and from 
this result the negative current could in principle be obtained. (Along the way, we also 
discovered a very concise derivation of the arrival time formula with a complex potential, 
Eq.dSTD). 

In Section 4 we looked at the particular case of measurements of position at two moments 
of time for initial states which are eigenstates of the flux operator. This includes in particular, 
states with negative eigenvalues which are clearly backflow states. We showed that certain 
probabilities have surprising non-classical properties and may be the basis of experimental 
test. We also noted that it is not just the backflow states with non-classical properties. 
Eigenstates with eigenvalues in the range (0, 1) can also be non-classical. 

In Section 5, we computed the current and the flux for states consisting of superpositions 
of gaussians. These states are important since they are experimentally realizable. Backflow 
is easily obtained with these states but the maximum amount of flux is very small, only 
about 16 percent of the maximum possible. 
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In Section 6, we looked for analytic expressions for states matching as closely as possible 
the numerically computed states giving maximal backflow, computed by Penz et al j4|. We 
first showed that the analytic expression conjectured by Penz et al, sinw^/w, has positive 
current at t = 0, contrary to the numerical solution, so cannot be the maximizing state. We 
then presented two candidate analytic expressions and computed the current at arbitrary 
times analytically. The plot of the current in each case had reasonably good agreement with 
the numerical solution, except at the end points t = ±1 of the backflow region. We computed 
the most negative flux of these states. In one case, the total flux is about 70 percent of the 
numerically computed maximum backflow, significantly better than any previous analytic 
expression for a backflow state. For this most negative flux state, we plotted the probability 
of remaining in a; < against time in Fig. (12). This gives a particularly striking illustration 
of the backflow phenomenon, showing a distinct period of increase in probability. 

In Section 7, we discussed the classical limit of backflow. The most interesting aspect of 
this is the issue, first noted by Bracken and Melloy, that the eigenvalues of the flux operator 
are independent of h. This appears to mean that there is a genuine quantum phenomenon, 
negative flux, which is independent of h and which does not appear to go away in the 
naive classical limit /I — )> 0. We showed that this situation starts to appear more physically 
sensible when the projectors used in the definition of the flux operator are replaced by 
quasiprojectors, which includes a physical parameter characterizing the imprecision of real 
measurements. The eigenvalues then do depend on h and the most negative eigenvalue 
becomes less negative as the quasiprojectors become more smeared. Furthermore, there is 
evidence that all the negative eigenvalues become zero or positive as /?. — > 0, restoring the 
naive classical limit. However, there are clearly more issues to explore around this question. 

The features of backflow elucidated here may be of value in designing experiments to test 
backflow. These and related ideas with be explored elsewhere. 
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